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Abstract 

Diffraction of electromagnetic plane waves by the gratings made by periodically corrugating 
the exposed planar boundaries of homogeneous, isotropic, linear dielectric-magnetic half-spaces is 
examined. The phase velocity vector in the diffracting material can be either co-parallel or anti- 
parallel to the time-averaged Poynting vector, thereby allowing for the material to be classified as 
of either the positive or the negative negative phase-velocity (PPV or NPV) type. Three methods 
used for analyzing dielectric gratings — the Rayleigh-hypothesis method, a perturbative approach, 
and the C formalism — are extended here to encompass NPV gratings by a careful consideration 
of field representation inside the refracting half-space. Corrugations of both symmetric as well as 
asymmetric shapes are studied, as also the diversity of grating response to the linear polarization 
states of the incident plane wave. The replacement of PPV grating by its NPV analog affects only 
nonspecular diffraction efficiencies when the corrugations are shallow, and the effect on specular 
diffraction efficiencies intensifies as the corrugations deepen. Whether the type of the refracting 
material is NPV or PPV is shown to affect surface wave propagation as well as resonant excitation 
of surface waves. 

PACS numbers: 42.25.Fx, 78.20.Ci 
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I. INTRODUCTION 



Diffraction gratings are not only exploited by nature for the production of color [l| but 
have also been extensively used in optics for several centuries j^j. During the last century, 
tremendous progress in manufacturing techniques made diffraction gratings as the spectral 
dispersive elements of choice. Spectacular progress was also made on the numerical solution 
of the time-harmonic Maxwell equations for diffraction gratings, thereby greatly facilitating 
analysis and design Given such remarkable developments, one would think that the last 
word on diffraction gratings is imminent; but the emergence of isotropic dielectric-magnetic 
materials exhibiting phase velocity vector opposed in direction to the time-averaged Poynting 
vector has opened new prospects for diffraction gratings. 

Typically, a diffraction grating is a slab of either a metal or a dielectric material whose 
exposed surface is periodically corrugated. When a plane wave is incident on this surface, 
it is reflected partially in the specular direction fixed by Snel's law jrll and partially in 
nonspecular directions fixed by the periodicity of the corrugated surface in relation to the 
free-space wavelength. In addition, specular as well as nonspecular refraction into the slab 
may also occur, depending on the type of material. Finally, reflections from the back surface 
also contribute to the overall reflection from the grating, but those need not be considered 
when examining the essence of the phenomenon of diffraction. As the angle of incidence is 
changed, specular as well as nonspecular components of the reflected field wax and wane, 
which phenomenon is technologically exploitable. 

What would happen if a diffraction grating were to be made of an isotropic dielectric- 
magnetic material? If the phase velocity and the time-averaged Poynting vectors in this 
material are co-parallel, then the effects are not qualitatively different from a diffraction 
grating made of simply an isotropic dielectric material while the quantitative differences 
are due to differences in the relative impedance and the wavenumber inside the diffracting 
material j^. This paper is devoted to the case when the phase velocity and the time- 
averaged Poynting vectors in the diffracting material are oppositely directed. Though several 
names have been proposed for this class of materials, we think that the most descriptive is: 
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negative phase- velocity (NPV) materials. In contrast, the phase velocity and the time- 
averaged Poynting vectors are co-parallel in positive phase- velocity (PPV) materials. PPV 
materials are, of course, commonplace. 

The plan of this paper is as follows: The boundary value problem of a diffraction grating 
is presented in Section II, with careful delineation of field characteristics in the refracting 
half-space. Three methods of solving the boundary value problem are extended in Section 
III to encompass NPV refracting materials. These methods are: (i) the Rayleigh-hypothesis 
method, (ii) perturbative approach, and (iii) the C formalism. Numerical results for corru- 
gations of both symmetric as well as asymmetric shapes are presented in Section IV, and 
the effects of replacing a PPV material by its NPV analog, or vice versa, are extracted from 
those results. An exp(— icut) time-dependence is implicit, with i = V— 1, u as the angular 
frequency, and t as the time. 



II. BOUNDARY VALUE PROBLEM 



In a rectangular coordinate system (x, y, z), we consider the periodically corrugated bound- 
ary y = g{x) = g{x + d) between vacuum and a homogeneous, isotropic, linear, passive, 
dielectric-magnetic material, with d being the corrugation period. The region y > g{x) is 
vacuous, whereas the material occupying the region y < g{x) is characterized by complex- 
valued scalars €2 = e2R + ie2i and /i2 = /^2/? + ^/^2/ that depend on a;. If this medium is of 



the NPV type, then the following three conditions hold equivalently 

(|e2| - e2R)(|/i2| - /i2ij) - e2//U2/ > 

e2RfJ'2i + e2/yU2ij < > • (1) 

e2i?|/U2| + /i2i?|e2| < 
None of the three conditions hold for a PPV material. 

A linearly polarized electromagnetic plane wave is incident on this boundary from the 
region y > g{x) at an angle 6^0, (|6'o| < vr/2), with respect to the y axis, as shown in Figure [T] 
This plane wave can be either s-polarized or p-polarized . Given the orientation of the 
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plane wave with respect to the grating plane (i.e., the xy plane), all reflected and transmitted 
fields must be linearly polarized in the same way as the incident plane wave [isj ]. 



y 




> X 



isotropic 
PPV/NPV 
material 



FIG. 1: Schematic of the boundary value problem. A plane wave is incident at an angle 
(l^ol < ^/2), with respect to the y axis on the periodically corrugated boundary y = g{x) = g{x + d) 
between free space and a homogeneous, isotropic, linear dielectric-magnetic material. 



(2) 



Let the function f{x,y) represent the 2;-directed component of the total electric field for 
the s-polarization case, and the z-directed component of the total magnetic field for the p- 
polarization case. After starting from the time-harmonic Maxwell equations, this function 
can be shown to be a solution of the Helmholtz equations 

(v2 + ^)/(a;,y) = 0, y>g{x, 

(v2 + #e2/^2)/(a;,2/) = 0, y < g{x, 

where c is the speed of light in free space (i.e., vacuum). Outside the corrugation region 
maxg{x) > y > mmg{x), the field f{x,y) is rigorously represented as a superposition of 
plane waves as follows: 

i{aox - (3^^^y) + 

i{anX + l3n\) , y>ma.xg{x) . (3) 
i {anX - (3^'' y) , y<ming{x) 



f{x,y) 



exp 



+ E 



oo (1) 
n=-ooCn' exp 



Eo< 
n 



oo (2) 

oo exp 



Here, Cn \ (— oo<n<oo;p=l, 2), are scalar coefficients to be determined by the solution 
of a boundary value problem, while 

a„ = ^ sin6'o + 2mx/d 
= 



- at 



> . 



-^€2 /i2 - a. 



(4) 



Both jdn'^ and (3n' are double-valued functions by virtues of their definitions as square-roots. 
If fin"* represents an upgoing plane wave, then — /^i'^'' represents a downgoing wave, and vice 
versa. We must choose the correct signs for all as upgoing plane waves as well as for all 
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f3n'^ as downgoing plane waves. If jSn' is real-valued, it corresponds to a propagating plane 



wave; otherwise, it indicates evanescence. 



Note that (3n^ is either purely real or purely imaginary; and the condition 



Re 


4"" 


,0] 








'0 




> Vn 


(5) 


Im 







is appropriate for upgoing plane waves in the vacuous half-space y > g{x). As the direction 
and/or the angular frequency of the incident plane wave change, (3n^ may change from purely 
real to purely imaginary, or vice versa. Such alterations are usually marked by noticeable 
discontinuities in the diffraction spectrums and, hence, are called Rayleigh-Wood anomalies 
although their occurences are not at all anomalous j^. 



The refracting half- 
by virtue of causality 



spa ce y < g{x) being filled by a material medium, 62/ > and 1121 > 



14l Iq . The refracted plane waves must attenuate as y ^ — 00, which 



requirement leads to the condition 



Im 



2) 



] > Vn . 



(6) 



This condition on Im 



(3n^ automatically fixes the sign of Re f3n' , regardless of the signs 



j(2) 



of e2R and ^2R', furthermore, the transformation {e2R — e2i?,/U2_R — /U2_r} alters the signs 



of the real parts of all /5n 



(2) 



In order to find the scalar coefficients Cn\ we must apply the boundary conditions at 
y = g{x). These conditions, expressing the continuity of the tangential components of the 
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total electric field and the total magnetic field, can be written as 

f{x,g{x)+) = f{x,g{x)-) 
fi ■ V/(x, g{x)+) = a-^n- Vf{x, g{x)-) 

where cr = ^2 for the s-polarization case and a = 62 for the p-polarization case, while n is 
a unit vector normal to the boundary y = g{x). 




III. METHODS OF SOLUTION 



Analytical solution of the stated boundary value problem is well nigh impossible, except 
in a perturbative sense. Early numerical techniques Q| relied on the so-called Rayleigh 
hypothesis, according to which the expansions © can be assumed valid at y = q(x)± [17] . 



Therefore, those techniques were not applicable for deeply corrugated boundaries 



limitations of the Ra yleig h hypothesis were overcome by the T-matrix method jsj, |l8| and 







m. The 



the C formalism [121 l2i|, of which the latter displays superior performance. All of these 
techniques require various degrees of computational proficiency, and so we applied three 
different methods of solution for NPV diffraction gratings — as checks on each other, as 
applicable. 



(a) Rayleigh-hypothesis method 

This method was enunciated by Rayleigh for gratings made by sinusoidally corrugating a 
perfectly reflecting sheet. According to his hypothesis, the expansions Q have to be used in 
the boundary conditions (|7j). Both resulting equations are then projected into the Rayleigh 
basis {exp(2 0:^3;)}^=^, in order to obtain a. linear system of equations for all Cn^ . The 

.(2) 



refraction coefficients Cn are then eliminated [21|, |22| to yield 
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(1) (2) ~ 2 



n=— oo 




ym 






(1-a) 


- 










Pm 








, where 












^mni^^) 


1 /"^ 

d Jo 


exp [- 


.27r 
-I —{m — 



D 



mO 



~Po ^ ~ Pm^ 



(1) _ nC^)] ^(1) 



X + ^^(^(a;)] . (9) 

" 

The summation on the left side of (jH)) has to be appropriately truncated, and the equations 
are then put in the form of a matrix equation which can be solved by standard numerical 



methods 



The Rayleigh hypothesis is valid when the corrugations are not deep, and t 



le limit of its 



applicability for sinusoidal gratings has been rigorously established. Millar 2J| showed that 
the Rayleigh hypothesis is applicable for perfectly reflecting gratings of sinusoidal shape with 



25|, who also 



maximum slope not exceeding 0.448. This limit was validated by Hill & Celli 
noted that the methods exploiting the Rayleigh hypothesis could yield acceptable results for 
steeper gratings. Depine & Gigli [26^ carried out detailed numerical studies to show that 
the Rayleigh hypothesis can be considered adequate for dielectric sinusoidal gratings with 
maximum slopes as high as ~ 0.92. 



(b) Perturbative approach 



A perturbative approach applies well when the corrugations are very shallow. The integrals 
Dmn{u) can be stated exactly as the power series 



Dmn (u) = '^ — u^g^^\m-n), 



where 



a 



2tt 

[g{x)Y exp(— im— x) dx 

(Jj 



(10) 



is the m-th Fourier coefficient of the function [(^(x)]-'. These coefficients can be obtained 
through the recurrence relation 



beginning with 



where Smn is the Kronecker delta. 



(12) 



(13) 



Assuming the expansion |22, |2 



-(1) = V 



j=0 



(14) 



we arrive at an iterative scheme, whereby the coefficient Cn'''\ j > 1, can be obtained in 



terms of all lower-order coefficients Cn''' ^\ ■ ■ ■ , Cn'^' as follows: 



.(i>o) 



q=l 



(15) 



This scheme commences with 



(1,0) ^ '^P" P" r 



(16) 



which is the planewave reflection coefficient for a perfectly flat boundary (i.e., g{x) = 0) 
and requires the computation of 

(/?J')/5P - aoan)(l - ^) + S(e2/X2 - a) 



and 



iP^J^^pi^^ + a^an)(l-a)-^(e2/X2-^) 



(17) 



Provided the series (jl4j) converges rapidly, the reflection coefficients Cn^ can be computed 
quite easily even on hand-held computers. 
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(c) C formalism 



In order to avoid the use of the Rayleigh expansions Q in the corrugation region, the C 
formahsm begins with the transformation 



V = y- g{x) . 

Accordingly, the Helmholtz equations change to 



-2g 



dxdv dv 



c 



where 



e{v) 



1 

£2 



> 
< 



and 



dx 



dx"^ 



(19) 



(20) 



(21) 



(22) 



Because the coefficients of the differential equation (jSHI) depend on f in a piecewise 
fashion, the f-dependence of / is of the form exp {ipv) in each of the two pieces f > 
and f < 0. Following references |3l and we expressed the x-dependences of g, f and 
df /dv in terms of Fourier series, and obtained the following matrix equation: 



+ 



[A] 



-2 



[O] 



-1 



(23) 



In this equation, [A] and [B^^^] are diagonal matrixes formed by a„ and , respectively; 
the (m, n) element of the Toeplitz matrix Q is the (m — ?T,)-th Fourier coefficient of g; [O] 
and [I] are, respectively, the null and the identity matrixes; while [J-'] and JF are column 
vectors formed by the Fourier coefficients of / and —idf/dv, respectively. Clearly, p^^ is an 
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eigenvalue of the 2x2 block supermatrix on the left side of (|23|) : and the eigenvalue spectrum 
of this supermatrix has to be determined for the regions above {p = 1) and below (p = 2) 
the corrugated surface v = 0. 



For numerical solution, the infinite system in ()23|) must be truncated. If only terms are 
kept in each Fourier series, each block in the 2x2 block supermatrix is a x matrix, thus 
resulting in 2A^ eigenvalues for each value of p. In each region, all eigenvalues not satisfying 
the radiation condition at infinity should be discarded in the representation of the diffracted 



field. Accordingly, for p = 1, only those eigenvalues are acceptable for which eit 



20 



ler p I S 



29|. 



real-valued and positive or p is complex-valued with positive imaginary part 
Similarly, for p = 2, when the region f < is filled with a PPV material, acceptable values of 
p must be either real-valued and negative or complex-valued with negative imaginary part. 
However, when the region f < is filled with a NPV material, acceptable values of p must 
be either real- valued and positive or complex-valued with negative imaginary part. Actual 
diffracting materials must be dissipative; hence, whether the region w < is occupied by a 
PPV or a NPV material, the criterion 

Im [p] < (24) 

suffices for p = 2. This criterion for the eigenvalues, together with the criterion (jH)) for 
selecting /?„ , are the central modifications that we have incorporated in the conventional 
C formalism for making it applicable to diffraction by either PPV or NPV corrugated half- 
spaces. 

For PPV materials, Chandezon et al. j^l have shown numerically and Li jsol has shown 
analytically that the real-valued eigenvalues and the lower-order complex-valued eigenvalues 

mm— mm. f'9^ 

of ()23j) converge to ±B^' as A^ increases. This property must also hold for NPV materials, 
since, as noted by Li [20|, a plane wave is an eigensolution of ^ and the transformation (fTT^ 
does not change the relevant eigenvalue, whether the refracting material is of the PPV or the 
NPV type. Indeed, indicates that when the signs of both and p2R are changed, the 
eigenvalues in the refracting half-space transform into their own complex conjugates. This 
is demonstrated by the sample results presented in Table 1. 
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cpluj 

r 1 


cp/uj 




n 


with = 11 


with = 25 






2.43796^20.14356 


2.43796^^0.14356 


2.43796 ±20.14356 





-2.43796 ±iO. 14356 


-2.43796 ± zO. 14356 


-2.43796 ±z0.14356 





2.11442^^0.16492 


2.10547^^0.16623 


2.10547 ±20.16623 


2 


-2.11442 ±i0.16492 


-2.10547 ±20.16623 


-2.10547 ±20.16623 


2 


1.71339 =F«0.20374 


1.71413 =F 20.20419 


1.71413 ±20.20419 


3 


-1.71339 ±i0.20374 


-1.71413 ±i0.20419 


-1.71413 ±zO. 20419 


3 


0.94258 T zO.30709 


1.00455^20.34841 


1.00455 ±20.34841 


4 


-0.94258 ± i0.30709 


-1.00455 ±i0.34841 


-1.00455 ±i0.34841 


4 


-1.00099 ±i0.34382 


-1.03946 ±20.33671 


-1.03946 ±i0.33671 


-5 


1.00099 TzO.34382 


1.03946^20.33671 


1.03946 ±i0.33671 


-5 


1.71661 TzO.19773 


1.73180^20.20210 


1.73180 ±i0.20210 


-4 


-1.71661 ±i0.19773 


-1.73180 ±20.20210 


-1.73180 ±z0.20210 


-4 



Table 1: Some eigenvalues of (j23p for a dielectric-magnetic material 
(e2 = ±6 ±20.1, p2 = ±1 ±20.1) computed for truncation parameters = 11 and = 25, 

in comparison with j3n , for h/d = 0.1, ood/c = 27r/0.5 and = 15°. The first three 
columns span upgoing and downgoing waves, but only downgoing waves are acceptable in 
the present situation. Therefore, acceptable values of and p must conform to the 

restrictions © and ((211), respectively. 



Once the foregoing changes have been incorporated, implementation of the C formalism 
proceeds as usual. As this is well-documented in the literature jigl . 20, [2^, we only sketch 
the procedure here for completeness. The following two steps are undertaken: 
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(i) First, the function f{x,y) is written as 

f f >g{x) 

fix,y)={ ^ y{ ■ (25) 

[<g{x) 

Here, the sectional field functions 

= exp[z{aoX-(3^'^y)]+ c« exp [z + + 

5^exp(za„x) ^ C«/«exp{zp«[y-^^(x)]} (26) 



and 



f('\x,y) = ci^^ exp Kx - /?f y)] + 

^exp(za^a;) ^ C^f exp {^pf [|/ - ^?(x)]} , (27) 

contain Cg^^ as unknown scalars with the index q indicating the g-th eigenvalue 
of the 2x2 block supermatrix in (j23|) . and fml are the successive elements of the top 
half of the corresponding eigenvector. The set U^*') contains the indexes corresponding 
to physically acceptable propagating plane waves, and we note that the set U'-^-* is 
always an empty set when the refracting material is dissipative. In contrast, the set 
y^P^ contains indexes corresponding to physically acceptable evanescent plane waves. 

(ii) Second, ()2(j|l and ()27|1 are rewritten in terms of the variables x and v only and then 
introduced in the boundary conditions ((Zj). A complete set of linear algebraic equations 
is thereby obtained for the sets of the unknown scalars Cn^ and c!f\ The 2N scalars 
are then calculated using standard methods 23^ . 

The C formalism, not invoking the Rayleigh hypothesis and therefore not limited to grat- 
ings with shallow corrugations, is a very efficient and versatile theoretical tool for modeling 
the electromagnetic responses of diffraction gratings of arbitrary permittivity and corrugation 
shape. As stated in Ref. [2^, the most distinctive feature of this formalism is its virtually 
uniform convergence, regardless of the incident polarization state and the permittivity of the 
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refracting material. Originally set up for uncoated, perfectly conducting gratings in classical 
mounts , the essence of the formalism — mainly, the simplicity of the coordinate system 
(fTI?|) — has allowed its extension to many other situations. Examples include multila yer- 
coated dielectric and metallic gratings \29i . conical mounting s 131 1. nonlinear materials 

nnn ft 

anisotropic materials [221 u3, 12^ nonhomogeneous materials j2a| , and crossed gratings 



32|, 
37|. 



This versatility of the C formalism is not matched by any other rigorous method for gratings 
[sj. We found that the characteristic features of this formalism are valid even for diffraction 
gratings of isotropic NPV materials. 

(d) Conservation of energy 



Diffraction efficiencies 



Re 



ci^f, (28) 



are defined for the propagating planewave components of the refiected field in the region 
y > maxg{x).. The normalized power Pa transferred across one period of the corrugated 
boundary into the refracting half-space y < g{x) can be calculated by virtue of the Poynting 
theorem, if the fields at the surface y = g{x)— are known. The principle of conservation of 
energy requires that 

Pa+ Yl < = (29) 

with Pa being completely absorbed by the refracting material. When we implemented any 
of the methods of solution presented in Sections III(a)-(c), we checked that the condition 
()29|) was satisfied to an error of 10 ppm. 

IV. RESULTS AND CONCLUSIONS 

Although corrugations of different shapes are used, we confined ourselves chiefiy to the 
most most widely used shape: 

g{x) = h cos \ —r-j ■ (30) 



d 
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Calculations of the diffraction efficiencies were made for many values of the geometric ratio 
h/d and normalized periodicity uod/c, using one or all three of the methods of solution 
described in Section III, as applicable. These results are presented and discussed in Sections 
IV(a)-(c). Asymmetric counterparts 

5f(x) = /ii cos(^^) + /i2 cos(^^ - 7) (31) 

of the symmetric gratings (plj) are addressed in Section IV(d). 



(a) Shallow gratings 



Let us begin with gratings described by (pH) . When the boundary y = g{x) is perfectly flat, 
the only non-zero reflection coefficient is The transformation 

{e2 -e*, H2 -/i2} ' (32) 

which amounts to the replacement of a NPV/PPV half-space by an analogous PPV/NPV 
half-space, changes the phase of Cq^"* but not its magnitude |l^; hence, the transformation 
does not affect Cq at all. 

For a shallow grating, we therefore expect that the magnitude of the specular reflection 
coefficient would not be greatly affected by the transformation ((221), but the effect of the 
transformation should be unambiguously evidenced by the nonspecular diffracted orders. 
This is indeed true, as was borne out by results computed using the perturbative approach 
of Section Ill(b). 

Figure El presents the diffraction efficiencies Cq and e^^ as functions of 9o G {— 11/2,71/2) 
when h/d = 0.07 and ud/c = 27r/0.8. The refracting material is of either the PPV (62 = 
5 + iO.Ol, /i2 = 1 + «0.01) or the NPV (e2 = -5 + iO.Ol, fi^ = -I + iOM) type. Calculations 
were made for both the s- and the j9-polarization cases. Two Rayleigh-Wood anomalies 
occur at ^0 ~ 11-54° ^ = 0) and at Oq ^ 36.87° = 0). 

Clearly, Figure |21 shows that the transformation ()32|) does not greatly affect Cq, except at 
low l^^ol- In contrast, the same figure shows that the nonspecular diffraction efficiency e^_^, 
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which is non-zero only for sin 6*9 > —0.2, is gravely affected by the type of the refracting 
material. 

The diversity can be understood as follows: When the boundary is perfectly flat, the 
transformation (jH^ leaves the magnitude of the reflection coefficient unchanged only for non- 
evanescent incident plane waves; but that is not a true statement for incident evanescent 
plane wave. Q. In the trough, of a shallow ,at,„g, ,.e., fo. ,na.„.) > , > the 
total field actually has both specular {n = 0) and nonspecular {n ^ 0) components, by 
virtue of the Rayleigh hypothesis. The nonspecular components are like evanescent plane 



waves because they are characterized by Re 



(1) 

'n 



0. Their presence ensures that the 



nonspecular reflection efficiencies, although weak for very shallow gratings, are considerably 
affected — in contrast to the specular reflection efficiency — by the transformation of the 
refracting material from the NPV/PPV to the PPV/NPV type. 
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FIG. 2: Diffraction efficiencies (a) Cq and (b) e^_i as functions of the incidence angle 9q, when 
h/d = 0.07 and ujd/c = 27r/0.8. The refracting material is of either the PPV (e2 = 5 + zO.Ol, fi2 = 
1 + iO.Ol) or the NPV (62 = — 5 + iO.Ol, /X2 = — 1 + iO.Ol) type. Calculations were made for both 
the s- and the p-polarization cases. Note that 69(^0) = Co(— ^0) and that e^_ii{6Q) 7^ only for a 
Umited 0o~range. The same results were obtained with all three methods of solution described in 
Sections III(a)-(c). 



(b) Deep gratings 
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As the corrugations grow deeper (i.e., as h/d increases in value), the transformation of the 
refracting medium from NPV/PPV to PPV/NPV increasingly affects the specular efficiency 
Cq as well. This is demonstrated by the plots of the diffraction efficiencies Cq _^ _2 versus 9q 
in Figure 01 for h/d = 1. As the Rayleigh hypothesis is inadequate for sinusoidal gratings 
with h/d^ 0.3 34], the presented plots were obtained using the C formalism. Incidentally, 
Rayleigh-Wood anomalies are evident in this figure at ^0 = 0° = 0) and at = 30° 
(/3«=/3« = 0). 




15 30 4S 60 75 90 1 5 30 45 60 75 90 

angle o1 incidence (deg) angle o1 incidence (deg) 



FIG. 3: Diffraction efficiencies (a) Cg, (b) e1^, and (c) e'^_2, and (d) normalized absorbed power Pa, 
as functions of the incidence angle when h/d = 1 and ujd/c = 27r/0.5. The refracting material 
is of either the PPV (e2 = 6 + iO.Ol, /xa = 1 + iO.Ol) or the NPV (e2 = -6 + iO.Ol, /ig = -1 + ^0.01) 
type. Calculations were made for both the s~ and the p-polarization cases, using the C formaUsm 
with N = 29. Note that e5((9o) = el^i-do), and that eLi{0o) / as well as eL2{0o) + only for 
limited ^o^ranges. 
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(c) Resonant surface-wave excitation 



A mechanism that can introduce dramatic changes in the diffraction efficiencies when the 
type of the refracting material is transformed from NPV/PPV to PPV/NPV is the resonant 
excitation of surface waves. Surface waves are not allowed to propagate on a plane boundary 
between vacuum and a material whose permittivity and permeability have positive real parts 
(Boardman 1982). For p-polarized (resp. s-polarized) surface waves to propagate along that 
boundary, the real part of the permittivity (resp. permeability) of the refracting material 
must be negative. Dielectric materials with negative real permittivity are exemplified by 
plasmas as well as metals (Boardman 1982). With the emergence of NPV materials, the 

jropagation of both types of surface waves on the same plane boundary has become possible 

2^|40|, although in principle for different frequencies. 

If the electromagnetic fields of the surface wave on each side of a plane boundary are 
described by (jHl) without the term corresponding to the incident plane wave and with only 
the n = terms in the two series, dispersion relations can be easily obtained js^ . Thus, 
the wavenumber ao of the surface wave satisfies the relation 

"° ^ /xi - 1 ' ^ ' 

for s polarization, and 

2 

"0 = ^^137^2, (34) 
for p polarization. These relations apply rigorously to plane boundaries only. 

To illustrate how the surface wave mechanism can affect the diffraction efficiencies of 
a grating when the type of the refracting material is transformed from NPV/PPV to 
PPV/NPV, even for shallow corrugations, we performed calculations with €2 = —1.8 + i0.01 
and fi2 = 1-5 + iO.Ol. According to the conditions (|H), this material is of the NPV type. 

We see, from (jSH), that a plane boundary can support a p-polarized surface wave with 
cao/uj ^ 1.63 + iO.Ol. If the transformation ()32|) is implemented, the plane boundary can 
not support the propagation of a p-polarized surface wave; instead, as follows from (j33|l . an 
s-polarized surface wave can then propagate with cao/u ~ 1.99 + i0.02. As the real parts of 
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both values of cao/uj are greater than unity, neither of the two types of surfaces waves can 
be resonantly excited by illuminating the plane boundary by a plane wave from the vacuum 
side. 
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42 1, surface waves can be 



However, as is well-established in the grating literature 
coupled to propagating waves through the periodicity of a corrugated boundary. If the 
period of the grating is convenently chosen, the surface wave can be coupled with one of 
the nonspecular components (i.e., n ^ 0). For example, after choosing ud/c = 27r/1.51 
and assuming that the wavenumber of the surface wave is not appreciably altered by the 
corrugation, (jU for a„ predicts a coupling at 6'o ~ 7°, when the refracting material is of the 
NPV type (e2 = —1.8 + zO.Ol, /i2 = 1-5 + zO.Ol). This is confirmed by the numerical results 
shown in Figure^ for h/d = 0.07. The zeroth-order efficiency curve as a function of angle of 
incidence (Figure Et^) for s-polarization is almost flat, whereas for p-polarization it exhibits 
a pronounced dip, near an angle of incidence very close to that predicted by the quasiplane 
approximation (jMj) for surface-wave excitation. This dip, at 6o ~ 7.9°, is not related to a 
redistribution of the incident power between other reflected orders (i.e., a Rayleigh-Wood 
anomaly). Instead, this dip is associated with a peak in the power absorbed in the refracting 
material, as can be seen from the Pa^Oo curve in Figure Eb- At the dip, nearly 87% of the 
]9-polarized incident power is absorbed by the refracting material, whereas just less than 2% 
of incident power is absorbed at all angles of incidence for the other polarization case. 



That the transformation of the type of the refracting material from NPV to PPV radically 
alters the conditions for surface-wave excitation is evident on comparing Figures ^ (e2 = 
1.8+iO.Ol, fi2 = -1.5+iO.Ol) andEl(e2 = -1.8+iO.Ol, fi2 = 1.5+iO.Ol). Three differences are 
noticeable. First, the polarization-dependences are different: whereas the the NPV grating 
exhibits a strong absorption peak for p- but not for s-polarization, its PPV analog exhibits 
a strong absorption peak for s- but not for p-polarization. Second, the absorption peaks 
occur at different angles of incidence for the two types of materials. The peak absorption 
in Figure 1^ occurs for p-polarization at Oq ~ 7.9°, but for s-polarization at 6o ~ 30.9° in 
Figure Eb- Third, although both peak absorptions are very strong, that in Figure near 
6^0 ~ 30.9°, is almost total (nearly 99% of the s-polarized incident power). 
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FIG. 4: (a) Diffraction efficiency Cg and (b) normalized absorbed power Pa functions of the 
incidence angle Oq, when h/d = 0.07 and tod/c = 27r/1.51. The refracting material is of the NPV 
(e2 = —1.8 + iO.Ol, iJ,2 = 1.5 + iO.Ol) type. Calculations were made for both the s- and the 
p-polarization cases using all three methods presented in Section 111. 
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FIG. 5: Same as Figure 0J except that the refracting material is replaced by its PPV analog 
(£2 = 1.8 + iOm, fi2 = -1.5 + iO.Ol). 

For gratings with deep corrugations, the wavenumber of the surface wave can be apprecia- 
bly different from the values in the absence of the corrugations, or the surface wave can even 
be forbidden to propagate. This can be concluded from Figures E] and Figures [71 which were 
drawn for the same parameters as for Figures |3] and Figures except that now h/d = 1. 
Apparently, the p-polarized surface wave does still play a role in the diffraction by the NPV 
grating (FigureEfc), with a broad absorption peak at 9o ~ 10.5°, close to the value found for 
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h/d = 0.07 in Figure |3)d. But no resonant behavior can be seen in Figure [71 the refracting 
material then being of the PPV type. The Rayleigh-Wood anomaly determined by [3^1 = 
is indicated in both figures at ~ 30.5°. 



1.0 - 




15 30 ^ 60 75 90 

angle o1 incidence (deg) 



FIG. 6: Diffraction efficiencies (a) Cg and (b) e^^ and (c) normalized absorbed power Pa as 
functions of the incidence angle when h/d = 1 and ud/c = 27r/1.51. The refracting material is 
of the NPV (e2 = — I.8 + i0.0I, = 1.5 + iO.Ol) type. Calculations were made for both the s- and 
the p-polarization cases using the C method with N = 29. 

(d) Asymmetric corrugations 

In order to illustrate effect of the corrugation shape on the diffraction efficiencies, we also 
considered asymmetric corrugations described by (jHT|) . Diffraction efficiencies were calculated 
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FIG. 7: Same as Figure IHl except that the refracting material is replaced by its PPV analog 
(e2 = 1.8 + iO.Ol, fi2 = -1.5 + iO.Ol). 

for hi/d = 0.12, h2/d = 0.078 and 7 = n/2, so that {m8ixg{x) — mmg{x)) /d ^ 0.33. 
Calculated values of e^, e!l^, e'^_2, and Pa as functions of the incidence angle 9q are plotted 
in Figure IHl for ood/c = 27r/0.5, when the refracting material is of either the PPV (e2 = 
5 + zO.Ol, /i2 = 1 + iOM) or the NPV (ea = -5 + zO.Ol, /ig = -1 + iOM) type. These 
plots were made for both the s- and the p-polarization cases, using the C formalism with 
= 29. Clearly, application of the C formalism is not limited to simple sinusoidal gratings. 
Additionally, as in Sections IV (a) and (b), the differences between NPV and PPV gratings 
are easy to divine from Figure |H1 and Rayleigh-Wood anomalies are present therein. 

The corrugation shape should affect surface wave propagation. This conjecture was verified 
when the calcuations for Figures El (e2 = —1.8 + zO.Ol, fi2 = 1.5 + iO.Ol) and El (62 = 1.8 + 
iO.Ol, /i2 = — 1.5 + iO.Ol) were repeated, but for the shape delineated by (jSH). Figures 01 and 
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FIG. 8: Diffraction efficiencies (a) Cg, (b) e^^^, and (c) elg, and (d) normalized absorbed power 
Pa, as functions of tlie incidence angle Oq, when ud/c = 27r/0.5. The corrugation shape is given by 
(jnU, with hi/d = 0.12, /12/d = 0.078, and 7 = 7r/2. The refracting material is of either the PPV 
(£2 = 6 + iOm, /i2 = 1 + «0.01) or the NPV (e2 = -6 + iO.Ol, /xs = -1 + iO.Ol) type. Calculations 
were made for both the s- and the ^-polarization cases, using the C formalism with N = 29. 

ITUl were drawn for hi/d = 0.04, h2/d = 0.026 and 7 = 7r/2. As {m8ixg{x) — ming{x)) /d ~ 
0.11 is rather small, the wavenumber of the surface wave should be predicted reasonably well 
by (j33j) and (j34|) . But the introduction of a Fourier harmonic to a sinusoidal corrugation 
seems to change strongly the coupling between the surface and incident waves, both for NPV 
and for PPV materials, as can be observed from Figures M and HUl 

Finally, we must remark on a major difference and a major similarity between Figures 
I2HZI on the one hand and Figures IHHTUl on the other. The corrugation shape is symmetric 
for the former set of figures, but asymmetric for the latter. We see that Pa{Oo) = Pa{So) 
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FIG. 9: Diffraction efficiencies (a) Cq and (b )e^]^, and (c) normalized absorbed power Pa as 
functions of tlie incidence angle when ivd/c = 27r/1.51. The corrugation shape is given by 
(PTT|) . with hi/d = 0.04, /12/d = 0.026, and 7 = 7r/2. The refracting material is of the NPV 
(e2 = —1.8 + iO.Ol, iJ,2 = 1-5 + iO.Ol) type. Calculations were made for both the s- and the 
p-polarization cases, using the C formalism with N = 29. 

for symmetric corrugations, but Pa{6o) 7^ Pai—60) for asymmetric corrugations. However, 
whether or not the corrugations are symmetric, eo(6'o) = eo(— ^o); as expected from reci- 
procity arguments (Petit 1980). 



V. CONCLUDING REMARKS 



In the foregoing sections, we extended Rayleigh's own method, a perturbative approach, 
and the C formalism to encompass diffraction by surface-relief gratings made of made of 
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FIG. 10: Same as Figure IHl except that the refracting material is replaced by its PPV analog 
(£2 = 1.8 + iO.Ol, = -1.5 + iO.Ol). 

isotropic, negative phase-velocity materials. This was enabled by carefully representing the 
field inside the refracting material. Numerical results for corrugations of both symmetric 
as well as asymmetric shapes were obtained, as also for both s- and p-polarized incident 
plane waves. We concluded that replacement of a PPV grating by its NPV analog affects only 
nonspecular diffraction efficiencies when the corrugations are shallow. When the corrugations 
deepen, the specular diffraction efficiencies are also affected by the type of the refracting 
material. 

Surface wave propagation as well as the resonant excitation of surface waves also depends 
on whether the refracting material is of the NPV type or its PPV analog. Excitation of a 
surface wave through a grating plays a fundamental role when high selectivity is desired. In 
common PPV gratings, surface waves have been exploited for efficient conversion of p- to 
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s- polarizations, or vice versa in conical mountings as well as to obtain enhanced nonlinear 
optical effects through the enhancement of surface fields usually associated with the resonant 



excitation of surface polaritons 
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42j |. Surface waves play an important role in the concept 



of a perfect lens realized using a NPV material, since the field of an image, which can not 
be focused by a normal lens, can be transferred through a NPV layer by the excitation of 



surface waves at both of its boundaries l6 



40|. We expect that NPV slabs with periodically 



corrugated boundaries combine both attributes. 
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